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ROLE OF CIRCUMFERENTIAL FLOW IN THE STABILITY OF 

FLUID-HANDLING MACHINE ROTORS 

D.E. Bently and A. Muszynska 
Bently Rotor Dynamics Research Corporation 
Minden, Nevada 89423, U.S.A. 

The recent studies of the dynamic stiffness properties of fluid lubricated bearings 
and seals by the authors [Refs. 1-7] have yielded most of the generalized charac- 
teristics discussed and used in this paper. They include bearing and seal nonlinear 
fluid film properties associated with stiffness, damping, and fluid average circum- 
ferential velocity ratio. Analytical relationships yield the rotor system's 
dynamic stiffness characteristics. This paper shows the combination of these data 
to provide the fluid- induced rotor stability equations. 
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Precession frequency at threshold of stability 
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RESULTS OF PERTURBATION TESTING OF ROTORS IN BEARINGS AND SEALS 


Dynamic stiffness properties of rotor systems are studied by applying a perturbation 
force of sinusoidal nature to a rotor. The best results are obtained when the force 
is circular (sinusoidal in two orthogonal lateral directions) and applied in a plane 
perpendicular to the shaft axis. The rotor is held at a constant speed while the 
perturbation force frequency is varied from zero to frequencies well above the 
rotative speed, in either forward (the same direction as shaft rotation) or reverse 
direction (opposing rotation). The frequency, amplitude, and phase of the perturba- 
tion force is accurately documented. The amplitude and phase of the displacement 
response of the rotor at each perturbation force frequency is then carefully mea- 
sured. Dividing the perturbation force vectors by the corresponding response 
vectors yields the dynamic stiffness vectors. Plotting these dynamic stiffness 
vectors in terms of their direct and quadrature components across wide perturbation 
frequency ranges allows dissection and evaluation of the dynamic stiffness compo- 
nents [Refs. 1,2, 4, 5]. 

In classical rotor dynamics for laterally isotropic systems, the bearing or seal 
fluid film direct dynamic stiffness terms are referred to as combination of direct 
(radial) spring, cross damping, and direct mass. The quadrature dynamic stiffness 
contains cross spring, direct (radial) damping, and cross mass. (Laboratory obser- 
vations up to this time, however, reveal no cross mass term.) 

The authors use the same definitions, but there are the following adjustments: (a) 
direct spring is shown to be a sum of the fluid film hydrodynamic stiffness plus the 
hydrostatic (externally pressurized) stiffness minus a component associated with the 
fluidic inertia effect; (b) cross damping is shown to be a component of the fluidic 
inertia effect; (c) the direct mass term is a component of the fluidic inertia 
effect; and (d) cross spring is a product of the direct damping, circumferential 
average velocity ratio, and rotative speed. 

For a particular seal or bearing operated at low eccentricity at a particular 
rotative speed, with a specific fluid at a known temperature (for knowledge of 
viscosity), the results of the combination of a forward and reverse perturbation 
runs typically yield the direct dynamic stiffness and quadrature dynamic stiffness 
curves as shown in figure 1. Additionally, a steady-state load test can be used to 
better evaluate the stiffnesses at a zero perturbation (precession) frequency (uj=0) 
[Ref. 2]. 


The general expressions of the fluid dynamic force components applied to the rotor 
with a coordinate system of +Y up and +X right and counterclockwise (X toward Y) 
rotation are 
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For sinusoidal isotropic (laterally symmetric) response, due to circular perturba- 
tion force F pert 2 at frequency u>, the force balance relationship for a concen- 

tically rotating shaft at the bearing or seal is as follows: 
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F pert z = [V (w '^ )2M fl +jD(u) ' Afi)+K ' Mu,2+jD s u>]z (1) 

where: 

z=x+jy is rotor lateral response vector at frequency w. 

Fpert z P er ^ ur ^ a ^^ on f° rce vector of the sinusoidal circular nature. 

K„ is the sum of the hydrodynamic plus the hydrostatic (externally pres- 

surized) bearing/seal fluid film direct (radial) stiffness (the non- 
linear character of this term as a function of shaft eccentricity is a 
vital relationship for the rotor stability solution. 




FIGURE 1 TYPICAL SEAL OR BEARING FLUID FILM DYNAMIC STIFFNESS PLOTS VERSUS PERTUR- 
BATION FREQUENCY FOR FORWARD AND REVERSE PERTURBATION. DIRECT DYNAMIC 
STIFFNESS CONTAINS STATIC RADIAL (DIRECT) STIFFNESS AND INERTIA TERMS. 
QUADRATURE DYNAMIC STIFFNESS CONTAINS RADIAL DAMPING AND "CROSS SPRING" 
TERMS (SHAFT DAMPING NEGLECTED). 
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j i s V"l- 

uj is sinusoidal perturbation frequency (+ for forward, - for reverse). 

Later, this same term means also rotor precession rate when instability 
(stable limit cycle, fluid- induced whirl or whip) occurs. 

D is fluid film direct (radial) damping. 

\ is the average circumferential component of fluid flow ratio to rotative 

speed in the seal, bearing, or rotor periphery. It is also the rate at 
which the fluid dynamic force rotates. It has been described in prior 
publications [Refs. 5-11]. Most authors assume that this rate is 1/2 
[Refs. 14-16]. The major change and improvement documented by the au- 
thors is that A is not constant and, in particular, is a nonlinear, de- 
creasing function of the shaft eccentricity. 

Q is rotor rotative speed. 

M-, is a synthetic expression with units of mass, called the fluidic inertia 

effect. 

C is bearing or seal radial clearance, 

e is bearing or seal eccentricity ratio e=|z|/C. 

K,M,D are rotor first lateral mode modal (generalized) stiffness, mass, and 
and external damping respectively. 

From (1) the rotor system direct dynamic stiffness is as follows: 

K D = (F pert z /|z|) (cos a) = K b" <“>-AQ) 2 M fl + K - Mu) 2 
The quadrature dynamic stiffness is 

k q = -(F pert Z /Iz|) (sin a) = D(urAfl) + D g iu 

where a is angle between the perturbation force and rotor response vectors. 

Expanding the binominal terms of the fluidic inertia effect, it may be observed that 
the fluid- related dynamic stiffness elements are as follows: 

(a) K xx = K yy = K B ~ i » so that the possibility of a "negative spring" 

exists, as is noted by the authors in prior reports [Refs. 1-5], and is well 
known to designers and researchers of squeeze film dampers. 

(b) D = D = +2AflwM f , . It is important to note that this term is positive in 

aj yx I I 

nature, as is being now reported by all researchers. This has been reported as a 
negative term in the 1950 to 1970 era. 

(c) *^y = = The ,,cross " spring is the stiffness of the fluid wedge 

support term which, as can be observed, is the product of the fluid average cir- 
cumferential velocity Afi, times the fluid direct damping D. It results from the 
rotative character of the fluid dynamic force [Refs. 5,11,14-16]. The term AftD 
is the Key of the rotor stability algorithm. 

Two vital characteristics of the fluid inertia effect are 

(1) The value of the fluidic inertia effect is zero when the precession speed 
u>=Afl. As most everyone knows, this corresponds, or is near, to one of the condi- 
tions of forward fluid- induced instability such as whirl. At the threshold of 
whirl, the fluidic inertia effect term is, therefore, quite ineffective. 
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(2) The fluidic inertia effect has the "ghost-like" characteristic because it 
disappears with increased shaft eccentricity following increasing radial load. 
Specifically, it is active only as long as circumferential flow is strong and the 
force from the "cross spring" term is dominant. This occurs at low through 
medium shaft eccentricities. At high eccentricities the circumferential flow is 
suppressed, yielding priority to secondary flows (axial and backward [Ref. 10]). 
At shaft high eccentricities the direct dynamic stiffness terms, mainly the 
hydrodynamic spring, become dominant and the fluidic inertia effect vanishes. 

Due to these characteristics, the fluidic inertia has very little effect on stabil- 
ity threshold or the orbit amplitude and precessional frequency of whirl or whip 
limit cycles. Therefore, it will be neglected in the following analysis. 


NONLINEAR CHARACTERISTICS OF BEARINGS AND SEALS 

The next required relationships for determination of rotor stability are nonlinear 
expressions for the fluid film direct stiffness Kg versus shaft eccentricity, direct 
damping D versus eccentricity, and the fluid average circumferential flow velocity 
ratio A versus eccentricity. It will soon be shown that the only vital nonlinearity 
affecting the size of whirl or whip limit cycle orbit (rotating eccentricity) is the 
direct stiffness term. Typical relationships are illustrated in figures 2, 3, and 
4. The nonlinear functions of the eccentricity ratio which are taken into consid- 
eration are as follows: 




(i- e 2 ) 3 


, K Bo = 4000 lb/in 


( 2 ) 



Ecc (e) 



FIGURE 2 TYPICAL FLUID FILM DIRECT 

(RADIAL) STIFFNESS VERSUS SHAFT 
ECCENTRICITY RATIO. 
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FIGURE 3 TYPICAL FLUID FILM DIRECT 
(RADIAL) DAMPING VERSUS 
SHAFT ECCENTRICITY RATIO. 



FIGURE 4 TYPICAL FLUID CIRCUMFERENTIAL 
AVERAGE VELOCITY RATIO VERSUS 
SHAFT ECCENTRICITY RATIO. 



FIGURE 5 ROTOR/BEARING MODEL. THE 
SIMPLEST SYSTEM DEMONSTRAT- 
ING WHIRL AND WHIP PHENOMENA. 


D 

D = > D q = 40 lb sec/in (3) 


A = A q (1 - e 2 ) 1 / 5 


(4) 


These functions qua! itatively exemplify the bearing/seal parameter behavior. They 
are not locked into a specific bearing/seal geometry, dimensions, or Sommerfeld 
number. Kg^, D 0 , and A 0 are fluid film radial stiffness damping and circumferential 

velocity ratio at concentric shaft respectively. 

The experimentally identified variation of A with eccentricity ratio is noted in the 
prior reports [Refs. 7, 9, 11]. No clean laboratory test has been found for obser- 
vation of A at this time, but a few indirect measurements have been accomplished. 
More numerical/analytical calculations, such as in [Ref. 10], will also help to 
identify the values of A as function of eccentricity and other parameters. 

ROTOR/BEARING SYSTEM 

The final set of equations for determination of stability is the equation of the 
rotor carrying bearing, seal, or other elements which provide the fluid. A simple 
rotor supported in one rigid and one fluid lubricated bearing is illustrated in 
figure 5. 

The dynamic stiffness, K , observed at the mass M of the rotor system in figure 5 
is eq 


K eq = Ki - u) 2 M + ju)D c + f 


K d K 2 


k b +k 2 


where Kg is fluid film complex stiffness: 


( 5 ) 
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( 6 ) 


K 8 = K b + jD(urAQ) 
Thus: 


K 2 (K 1 -u) 2 M)+K d (K 1 +K 2 -u) 2 M)-D D(u)-Afl)+j[u)D fK 2 +K D )+D(u)-Afi)(K 1 +K 2 -u) 2 M)] 

1/ _ D S S D /-y\ 

N eq " K 2 + Kg + jD(urAO) 

In Eq. (7) K X ,K 2 are shaft partial modal stiffnesses (first bending mode); M is 
rotor modal mass; D is shaft modal external damping. Fluidic inertia effect is ne- 
glected. More forirral analysis of this rotor/bearing system is in [Refs. 3 and 13]. 


ROTOR STABILITY ANALYSIS 

From the standpoint of stability, the numerator of K represents poles and the 

M 

denominator represents zeros. Only those poles are of interest for stability. 
Further, since only a steady level of the instability limit cycle orbit, if it 
exists, is of interest, the exponentially increasing and the exponentially decreas- 
ing orbits (transient processes) are not discussed. This allows the direct and 
quadrature parts of the numerator of dynamic stiffness K (the denominator of 

dynamic motions) to be separated and set to zero individually to find the stability 
roots. Another way to state this, is that the roots are evaluated at log decrement 
of zero. From control theory backgrounds, the root locus shown in figure 6 is 
solved only when the direct (real) parts of the roots are zero. The corresponding 
quadrature value is the precession rate, leading the system to the limit cycle of 
whirl or whip. 

Liapounoff contributed strongly to the study of stability. One of his theorems 
noted in [Ref. 17] is 

"Liapounoff 1 s theorem states that under certain conditions which are frequently 
encountered in physical problems, the information obtained from the linear 
equations of the first approximation is sufficient to give a correct answer to 
the question of stability of the nonlinear system." 

Restated, the linear portions of the direct and quadrature stiffness K gq solved for 

their roots yield the answer to the first problem of stability. Once the question 
of stability is decided, the next step is to evaluate the amplitude of the limit 
cycle orbit and rotor precession speed as a function of the nonlinear character- 
istics. 


The direct and quadrature components of Eq. (7) are equalized to zero and solved 
for Kg and Cl so that 


-K 2 (K 1 -u) 2 M)(K 1 +K 2 -u) 2 M)-u) 2 0 2 K 2 

(Ki+Kg-u^M) 2 + u) 2 0 2 


and 
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The right-hand side expression of Eq. (8) for various precession rates, omit 
ting the minor terms with D g , is illustrated in figure 7. 


EXAMPLES 

In order to plot the precession speed w against the rotative speed Q for a partic- 
ular case, it is necessary to choose values and functional relationships for the 
various rotor system parameters. The most simple process of solving the problem is 
to step the eccentricity ratio "e" gradually from 0 to about 0.9. For each eccen- 
tricity the values of Kg, D, and \ are calculated from appropriate equations such as 
Eqs. (2), (3), (4). Eq. (8) is then solved for the precession u> associated with 
corresponding eccentricity ratio e. Finally, for each rotative speed fi, the 
size of the limit cycle whirl /whip orbit eC (where C is bearing or seal radial 
clearance), and the speed of precession w, are all established. 

As an example, for CASE A the nonlinear terms are stated as follows: 

„ _ ^Bo , K„ = 4000 lb/in (fluid film radial stiffness (10) 

d ” (l-e^) 2 * B0 at concentric shaft) 



FIGURE 6 EXAMPLE OF ROOT LOCUS OF THE WHIRL/WHIP ROOT (u» . IS THRESHOLD OF STA- 
BILITY). 


I 
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D = D 0 = 40 1b sec/in (constant for this case) 
A = \ 0 = 0.48 (constant for this case) 


( 11 ) 

( 12 ) 


The remaining parameters are chosen as 

Ki = 2000 Ib/in M = 1.0 lb sec 2 /in 

(13) 

K 2 = 38000 lb/in D g = 4.0 lb sec/in 

The graph of the precession frequency of the whirl and whip against the rotative 
speed is shown in figure 8. This figure also shows the whirl and the whip ranges, 
and the smooth transition region between whirl and whip. It also shows a few 
corresponding whirl and whip orbits at the bearing (or seal) with their magnitude 
and frequency, as well as the stability threshold and asymptotic frequency. 



FIGURE 7 RIGHT-HAND SIDE EXPRESSION OF EQUATION (8) VERSUS PRECESSION RATE. 

EQUATIONS (10) TO (13) INDICATE VALUES OF PARAMETERS USED IN THIS 
EXAMPLE. 
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Figure 8 also illustrates the result of CASE B, where A is no longer considered 
constant, but A = 0.48 (1-e 2 ) 1 ^ 5 . The result is virtually the same as CASE A, with 
constant A. 

40 

For CASE C, damping is no longer constant but equals D = IT may be noted 

in figure 8 that the whirl /whip frequency of CASE C is much more like practically 
observed behavior, much closer to the asymptotic frequency, which is natural, 
because damping is always a nonlinear term. 

From Eq. (9) the ratio of the precession rate uj to rotative speed rate Q is an 
intrinsic function of u>: 



FIGURE 8 ROTATIVE SPEED VERSUS WHIRL/WHIP PRECESSION FREQUENCY FOR THE ROTOR 
ILLUSTRATED IN FIG. 1 (C = 5 MILS, SHAFT CENTERED). WHIRL/WHIP ORBITS. 
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(14) 



1 

1 + d((k 1 + k 2 -u>w + 


The shaft damping D has a fair amount of control over the actual precession fre- 
quency uj to rotative speed Q as shown in Eq. (14), but has very little control over 
the right-hand side expression in Eq. (8). Therefore, dropping shaft damping D 
from Equation (8) and solving the equation for precession rate w yields 


KiK 2 + K B (Ki+K 2 ) ** 

W * [_ (K b + K 2 )M J (15) 

The lowest precession rate u>^ that the whirl or whip can have, is at the minimum 
Kg, which is at Kg 0 for the case of the shaft centered in the seal or bearing con- 
sidered in this paper. Therefore, in terms of precession frequency w, the threshold 
of stability is 


u> 


th 


M 2 + k Bo (K!+k 2 ) 
(k Bo + K 2 )M 



(16) 


Poritsky [Ref. 18] nearly had this thres hold more than thirty years ago. He showed 
the rotative speed whirl threshold at 2 1 /K 1 /M. That is excellent for that time. Note 
that if bearing stiffness Kg = 0, shaft damping D g = 0, and X 0 = 1/2 and thus Cl = 

2w, it would be right on. His paper apparently was subjected to very heavy criti- 
cism, so his excellent work was not pursued. The authors privately call the zero 
bearing stiffness threshold the "Poritsky threshold." It is useful to observe that 
if the outboard shaft stiffness K x (the portion of shaft away from the fluid lubri- 
cated bearing) becomes very soft (approaches zero), whirling may occur at extremely 
low rotative speeds. The authors demonstrated a laboratory model with this phenom- 
enon for several years before fully understanding its rules of behavior. 

The asymptotic frequency u> ASY occurs as rotative speed increases so that the shaft 
orbit takes up most of there learance, and eccentricity ratio e approaches 1. The 
limit occurs not with nonlinear terms X or D, but with bearing stiffness Kg which 
goes to extremely high values at high eccentricity so that for Eq. (15) 


u> ASY = “I e-»l = V(Ki+K 2 )/M (17) 


It should be noted that shaft damping D was assumed to be 0 for Eq. (15). If 
shaft damping D is high enough, it is possible for the rotor system to escape 
whipping as shown in [Refs. 12 and 13]. Virtually most documented cases from 
actual machinery, however, exhibit the asymptotic behavior. 
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OBSERVATIONS AND CONCLUSIONS 


There are numerous straightforward observations available from these stability 
algorithms. 

First, while the "critical s" of machinery are often posted on the nameplate as a 
fixed frequency, it is obvious that they vary widely due to the nonlinear stiffness 
of fluid film in bearings and seals. The whirl/whip root of Eq. (8) shows clearly 
that the resonance ("critical") varies from the threshold to the asymptotic fre- 
quency u>^gY as a function of the nonlinear stiffness Kg. In the examples shown, it 

varies from the threshold of 75 rad/sec (716 rpm) to the asymptotic frequency of 
2 00 rad/sec (1909 rpm) (Figure 8). This fact will probably be met with considerable 
disbelief, until the nature of the fluid- induced whirl and whip instabilities become 
more widely understood. 

Since both of the quadrature dynamic stiffness elements contain the direct damping 
term: K Q = uuD - Af2D the magnitude of the direct damping D has very little effect 
on stability for rotors operating at low eccentricity. This means that large 
direct damping D does not create stability. 

Direct damping D is well known to be a function of bearing or seal geometry and 
physical dimensions. It also varies directly with lubricant viscosity q. However, 
it is not the viscosity in direct damping D which influences stability. It is the 
influence of viscosity in the hydrodynamic, as well as the hydrostatic (externally 
pressurized) portions of the bearing or seal stiffness term Kg that exhibits 
influence on stability. 

It is already known that anti swirling mechanisms which suppress the strength of 
circumferential flow and drive the fluid average circumferential velocity ratio X as 
nearly as possible to zero are the proper prevention of fluid-induced whirl and whip 
[Ref. 6]. 

Reasonably well known is that forward preswirling of the fluid will drive the pre- 
cession to rotation ratio to higher values [Ref. 6], sometimes well above half ro- 
tative speed. A report on a pump [Ref. 19] with a whirl which tracked about 80% of 
rotative speed was presented recently. Another verbal report [Ref. 20] on a double 
volute pump on an offshore platform noted a whirl which tracked 75% of rotative 
speed. A double driven generator connected by hydraulic clutches with two gas 
turbines exhibited subsynchronous vibrations with frequency 82.5% to 89.7% of 
the rotative speed when the clutches were not entirely engaged [Ref. 21]. An 

"inside-out" bearing with a hollow shaft rotating around a stationary post, will 
produce X of slightly over 1/2 [Ref. 22]. 

The pump with seal whirl [Ref. 19] was most unusual, because nearly all rotor/seal 
systems exhibit whip only. With good probability, this is because there is a much 
softer spring K t in that pump. Usually when the seal is involved, the stiffness K* 
is much higher. It may be noted in the discussed example that an increase of the 
magnitude of Kt leads to the increase of the inception threshold of stability. 

It is reasonably well known that any force which moves the rotor to high dynami c 
eccentricity, thus raising bearing or seal stiffness, can prevent whirl or whip. 
This action and the stability thresholds due to the unbalance force excitation were 
previously published [Ref. 3]. 
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It is also well known that steady-state sideload forces, which move the rotor to 
static high eccentricity in the seal or bearing, are a pallative for fluid- induced 
instability. Newkirk [Ref. 23] made the first pressure dam bearing for that purpose 
in 1934. Deliberate sideloading (often by gravity for horizontal rotors) is a very 
popular "fix" for stability. 

Two errors of first magnitude correlated to steady-state sideloads occur. The first 
is the rule that "the heavier the rotor, the better the stability," which has been 
published (and apparently believed) for over 30 years. As can be observed from the 
stability algorithms, the more the rotor mass, the worse the stability for the 
shaft-centered case. However, since the rule surely was meant to apply to hori- 
zontal rotors in the gravity of earth, having more mass means more weight, resulting 
in higher eccentricity, and thus higher Kg, yielding the better stability. (This 
assumes no other steady-state si deloading “forces, such as caused by misalignment, 
which adds or subtracts from the gravity load.) Two stabilizing and destabilizing 
factors, namely the weight-related preload and the added mass, are both very effec- 
tive, so that which factor wins depends upon the particular conditions in each 
particular case. A formula for stability that utilizes clearance and rotor weight 
should be carefully avoided, no matter how popular it has been. 

The second big problem of a steady-state sideload application occurs when the rotor 
is stacked (not integral), consisting of shrink-on wheels or, especially, when the 
rotor is bolted together (known as a shish ke-bab design). The bowing resulting 
from the steady sideloading strongly activates the hysteretic internal friction. 
Thus, attempts to sideload the rotor may help reduce the fluid-induced instability 
by raising K g , but get the rotor into an even worse problem of hysteretic friction 
instability, 1)1 us possible rotor-to-stator rubs. 

There exists another rule which is incorrectly applied by rotating machinery users. 
It states: "For stability maintain oil pressure low." This rule generally holds 
true when the bearing has a partial fluid film (180 degrees as opposed to 360 
degrees lubrication). As long as the partial film is maintained, a regular circum- 
ferential flow pattern is not established so that stability is assured. To maintain 
partial lubrication, the pressure should be kept low; otherwise the bearing becomes 
fully lubricated with stronger circumferential flow, which makes the rotor suscep- 
tible to instability. The rule, "Keep oil pressure low," became however, universal, 
and thus often incorrectly applied. Experimental evidences, as well as the con- 
clusions from the stability algorithm discussed in this paper, show clearly that 
fluid higher pressure , which directly increases the fluid film radial stiffness Kg, 
improves rotor stability. 

An interesting sidelight is that with the discussed stability algorithm it is simple 
to get some idea of how much torque is developed by the whirl /whip instability 
action. With a shaft radius R, and whirl or whip orbit radius eC: 

Wedge Stiffness AQD lb in Whirl/Whip Torque \Q0eCR lb in 

Tangential Force AflDeC lbs Whirl/Whip Power Loss AfiDeCRu) lb in/sec 

For A. = 0.5, Q = 250 rad/sec, D = 40 lb sec/in, e = 0.5, C = 0.01 inch, and R = 10 
inches, the whirl/whip torque is 250 lb/in, and the power loss at precession fre- 
quency 124 rad/s is 31000 lb in/sec. 

Finally, the fluid dynamic force relationships discussed in this paper can success- 
fully be used in modelling more complex rotor/bearing/seal systems [Refs. 9,12, 
13,22]. The model predicts existence of several whirl/whip regimes and explains 
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FIGURE 9 SPECTRUM CASCADE OF THE ROTOR VIBRATIONAL RESPONSE DURING START-UP, 
MEASURED BY THE PROXIMITY TRANSDUCER "B." THE PLOT SHOWS A "JUMP" OF 
SELF-EXCITED VIBRATIONS FROM WHIP MODE TO THE SECOND MODE WHIRL (JUMP IN 
AMPLITUDE AND FREQUENCY). THE SKETCH INDICATES TRANSDUCER LOCATIONS. 
OSCILLOSCOPE IN ORBITAL MODE FOR SIGNALS FROM TRANSDUCERS A AND B DISPLAYS 
THAT SHAFT AT LOW SPEED WHIRL VIBRATES IN PHASE. AT HIGH SPEED WHIRL 
SHAFT VIBRATES AT ITS SECOND MODE (RIGHT- AND LEFT-HAND SIDE SHAFT SEC- 
TIONS ARE 90 DEGREES OUT OF PHASE). 

the "jump" phenomenon observed on laboratory rigs and shown in figure 9, as well as 
in reported field data (figure 10). 

The role of the fluid circumferential flow in rotor stability is emphasized. The 
fluid average circumferential velocity ratio is the key factor of the "cross spring" 
term, which, in turn, represents a portion of the quadrature dynamic stiffness: Kg 

= D(urAfi). The value AA used to be referred to as "bearing resonance" across half a 
century. While this resonance is hidden in synchronous perturbation tests, it is 
clearly exposed by nonsynchronous perturbation testing [Ref. 1-5]. An increase of 
the circumferential flow strength leads to an increase of A, thus a decrease of the 
quadrature stiffness of the system. When at the same frequency both K Q and direct 
dynamic stiffness Kr. equal zero, the instability occurs. This was shown^on the root 
locus of figure 6. The quadrature dynamic stiffness, combined with the nonlinear 
bearing (or seal) stiffness in the rotor system equation, provides the reliable 
stability prediction formulas. 
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FIGURE 10 SPECTRUM CASCADE OF VIBRATION RESPONSE OF A STEAM TURBINE DRIVEN COMPRES- 
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